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We show that the physical subspace in the Z2-slave-spin theory is conserved under the time 
evolution of the system. Thus, when restricted to the physical subspace, this representation gives 
a complete and consistent description of the original problem. In addition, we review two known 
examples from the existing literature in which the projection onto the physical subspace can be 
relaxed: (i) the non-interacting limit in any dimension at half filling and (ii) the interacting model 
in the infinite dimensional limit at half filling. In both cases, physical observables are correctly 
obtained without explicit treatment of the constraints which define the physical subspace. In these 
examples, correct results are obtained, despite the fact that unphysical states enter the solution. 
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A defining property of so-called "slave-particle" con- 
structions is that the physical degree of freedom is rep- 
resented by auxiliary degrees of freedom in an enlarged 
Hilbcrt space. In the Z2-slave-spin representation intro- 
duced in Refs. 1 and 2, the local electron operator da is 
represented by an auxiliary fermion f ir7 and a spin- 1/2 
slave-spin U = (If , If, I?) as c ia = 2Iff ia . The physi- 
cal states |phys) form a subspace in the enlarged Hilbcrt 
space and are singled out by the local constraints 



is indeed the case: one can verify that 



A.lphys) = 0, Ai = I* + 1/2 - (m - 1)' 



(1) 



where n, = f\^fi\ + fj^fu denotes the local /-fermion den- 
sity and i = 1,...N S runs over all the lattice sites. In 
terms of these new variables, the original Hubbard model 
can be written as 

H S = -U £ (Ifl?flf jtr + h.c.) + ^£(1-2/?). (2) 

In Refs. 1 and 2 we argued that the original problem is 
equivalent to the slave-spin problem Eq. (2) acting on the 
physical Hilbert space defined by Eq. (1). Explicitly, the 
projection operator onto the physical subspace was given 
in Rcf. 2 as: 



v = Y[0-~Qi), Qi = A l 



(3) 



The operator Qi = Af fulfills [Hs,Qi] = in the ex- 
tended Hilbcrt space but does not provide an additional 
constraint on the physical Hilbert space. This follows im- 
mediately from the fact that Eq. (1) implies Qi|phys) = 0. 
Nevertheless, because Qi has eigenvalues and 1 it is 
used to construct the projector. 

In order that this projective scheme can be imple- 
mented consistently, the physical subspace needs to be 
invariant under the time evolution generated by H$- This 



B l := [A h H s ] 



8t£W(/^ + /,Ua) 



M (4) 



where the sum runs over all sites i connected to site I by 
the hopping amplitudes. Clearly, Bi + in the extended 
Hilbert space but for physical states, it follows that 



B ; |phys) = 0, 



(5) 



i.e., the physical subspace as defined by Eq. (1) is con- 
served under the evolution generated by H$- The slave- 
particle representation given by Eq. (2) acting on the 
physical Hilbert space therefore indeed gives a complete 
and consistent description of the original problem. This 
is the view-point we have adopted in our original work. 2 
A more complete discussion of constraint quantum sys- 
tems, as initiated by Dirac, can be found in the text- 
books, see e.g. Refs. 3 and 4. 

In approximative treatments of the slave-spin problem, 
the restriction to the physical subspace is no longer guar- 
anteed. In particular, if we replace the condition Eq. (1) 
by the weaker condition (Ai) = for states in the ex- 
tended Hilbcrt space, unphysical states are no longer ex- 
cluded from the theory. The objection formulated by the 
authors of Rcf. 5 is built on this observation but is more 
general: they seem to be concerned by the fact that the 
mean-field treatment does not correctly enforce the pro- 
jection onto the physical subspace. 

Here we would like to clarify a potential source of con- 
fusion. In our original article, 2 we adopted the projec- 
tive view-point of the slave-spin construction, as reviewed 
above. In Sec. VI we also discussed a mean-field treat- 
ment which enforces (Ai) = on average by the Lagrange 
multiplier method. We stress that this approximation 
(with a static Lagrange multiplier) was not meant as an 
alternative to the projective scheme - to obtain a phys- 
ical state, one would still need to project the mean-field 
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solution to the physical subspace. However, the inclu- 
sion of the Lagrange multiplier term allowed us to access 
magnetic phases on the mean-field level. 

Interestingly, the sole fact that a given slave-particle 
state is not an element of the physical subspace does 
not imply that expectation values of physical observables 
are incorrect. In the following, we review two known 
examples from the existing literature in which correct 
results are obtained even if the constraints are relaxed 
completely. The first example has already been given in 
our original work, 2 namely the non-interacting system at 
half-filling. For completeness, we review here some of the 
main arguments. Ignoring the constraint on the physical 
Hilbert space Eq. (1), a ground state of Eq. (2) with U = 
at half-filling is readily found: 

tf |$o) = £o|$o}, |$o) = |{l/2})/®|0 o )/- (6) 

Here, 7f|{l/2})/ = +l/2|{l/2}} 7 for all i and \<f> ) f is the 
half-filled Fermi sea obtained from Eq. (2) by replacing 
all operators If with 1/2 at U = 0. This ground state 
|$ ) fulfills 

(Ai)* o =0, Ai = If + 1/2- (m-lf, (7) 

on average; it is however not an element of the physical 
subspace because 

(Qi)*o=(A 2 i h = 1/2*0. (8) 

Yet, it can be easily shown 2 that all expectation values 
of physical observables are correctly obtained from |$ ). 
Hence, the correct physics in the non-interacting limit 
at half filling is obtained despite completely relaxing the 
constraint. However, if needed, one can obtain the phys- 
ical ground state by projection of |$o) and the result is 
a state which is a equal-weight superposition of all the 
states obtained from |$o) by a gauge transformation. 2 



Another example in which the projection onto the 
physical subspace can be ignored has been pointed out 
by Schiro and Fabrizio in Ref. 6 and elaborated on by 
Baruselli and Fabrizio in Ref. 7. These authors studied 
the single-impurity Anderson model as well as the half- 
filled Hubbard model in the infinite dimensional limit for 
which they derived the partition function in the slave- 
spin representation 

Z=(-J Z S , Z s = Tr(e-^). (9) 

Here, the interacting slave-spin Hamiltonian is given by 

Hs = ~^= E (lf/J/Uv + ^-) + TE(l-2^) J (10) 

and the coordination number z is taken to infinity at 
the end of the calculation. Remarkably, the projection 
to the physical subspace effectively drops out and enters 
the calculation of the partition function only through the 
prefactor (1/2) N (iV=number of sites) which accounts for 
the additional degrees of freedom in Hs- 

A general way to study if the projection to the physi- 
cal subspace qualitatively changes the physics associated 
with a given variational (mean-field) ansatz is to investi- 
gate the role of the "gauge fluctuations" . These are (un- 
physical) fluctuations of the mean-field parameters which 
enter the theory as a result of the gauge redundancy in- 
troduced by the slave-particle representation. For the Zi 
slave-spin theory, the role of the gauge fluctuations has 
been discussed in Ref. 8. It was also found that in the 
strongly correlated limit a fractionalized non- Fermi liquid 
phase with Z2 topological order and potentially gapless 
charge and spin excitations can exist. 
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